#include <math.h>
           
double j0(x);
double j1(x);
double jn(n,x);
double y0(x);
double y1(x);
            
double yn(n,x);
            
double x;           значение с плавающей точкой
int n;              целый порядок

Описание.
Процедуры j0, j1, jn вычисляют соответственно функции Бесселя первой степени порядка - 0, 1 и n.
Процедуры y0, y1, yn вычисляют соответственно функции Бесселя второй степени порядка - 0,1 и n. Аргумент x должен быть положительным.

Возвращаемое значение.
Эти функции возвращают результат функций Бесселя для x. Если x является отрицательным для y0, y1, yn, процедура устанавливает errno в EDOM, печатает сообщение об ошибке DOMAIN в stderr и возвращает отрицательное значение HUGE.
Обрaботкa ошибок может быть модифицировaнa при изменении процедуры matherr.

Пример:

            
#include <math.h>
            
double x, y, z;
.
.
.
y = j0(x);
z = yn(3,x);